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Abstract 



We present a comparison of numerical predictions on baryon asymmetry of the 
universe (Yb) for different neutrino mass models. We start with a very brief 
review on the main formalism of baryogenesis via leptogenesis through the de- 
^p , cay of heavy right-handed Majorana neutrinos, and then calculate the baryon 

f^ ' asymmetry of the universe for known six neutrino mass models viz., three quasi- 

■^ , degenerate, two inverted hierarchical and one normal hierarchical models, which 

are generated from canonical seesaw mechanism. The corresponding mass matri- 
ces for the right-handed Majorana neutrino Mjifi as well as the Dirac neutrino 
rriLR are fixed at the seesaw stage for generating correct light neutrino mass 
matrices niLL consistent with observations, and are again employed in the cal- 
culation of baryogenesis. This two-tier procedure removes possible ambiguity 



Qh! on the choices of m^ji and Man, and fixes the values of input parameters at the 



seesaw stage. We find that the ranges of predictions from both normal hierarchi- 
cal model (NHT3) and the degenerate model (DcgTlA) having mass eigenvalues 
(toi, — TO2, 7713), are almost consistent with the observed baryon asymmetry of 
k> I the universe. Other two degenerate models (DegTlB, DegTlC) and one inverted 

\^ • hierarchical model (InvT2A) lead to very small baryon asymmetry, Yb < 10~^^, 

5t ! whereas inverted hierarchical model (InvT2B) gives larger Yb > 10"^. Com- 

bining the present result with other predictions such as neutrino masses and 
mixings, and stability under radiative corrections in MSSM, the normal heirar- 
chical model appears to be the most favorable choice of nature. Possible sources 
of uncertainty in the estimations are pointed out. 
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1 Introduction 

Left-right symmetric GUT models such as SO(10)GUT, assume the existence of 
heavy right-handed Majorana neutrinos, and also predict small non-zero neu- 
trino masses through the celebrated seesaw mechanism[l]. Various neutrino 
mass models describing neutrino mass patterns, have been proposed[2] within 
the framework of seesaw mechanism, and their predictions on neutrino masses 
and mixings have been studied thoroughly. In our earlier communication [2] we 
had studied the stability criteria of these models in the presence of left-handed 
triplet Higgs field in type-II seesaw formula, and also the stability under ra- 
diative corrections in MSSM. The above criteria had been used to discriminate 
the neutrino mass models describing different patterns of neutrino masses. As a 
continuation of our earlier effort in this direction, we are now interested in the 
present work to have a comparison of the different theoretical predictions of the 
cosmological baryon asymmetry derived from the above neutrino mass models. 
This information may be treated as an additional point in the process of further 
discrimination of neutrino mass models. 

In SO(10)GUT, the role of heavy right-handed Majorana neutrinos is man- 
ifold. In addition to its role in seesaw mechanism, it also plays an important 
role to explain the observed baryon asymmetry of the universe [3], 

rl^Ms = (6.it!^i)io-i°. 

Such baryon asymmetry can be dynamically generated if the particle interac- 
tion rate and the expansion rate of the universe satisfy Sakharov's three fa- 
mous conditions [4]: (i) baryon non-conservation, (ii) C and CP violation and 
(iii) departure from thermal equilibrium. Being Majorana particles the heavy 
right-handed neutrinos automatically satisfy Sakharov's one of the above three 
conditions i.e., C and CP violation as they can have asymmetric decay to lepton 
and Higgs particles, and the processes occur at different rates for particles and 
antiparticles. Baryogenesis through leptogenesis by the electroweak sphaleron 
process has been widely accepted as a correct approach and a lot of work has 
been done in this direction[5,6,7,8,9,10]. The electroweak sphaleron process 
violates baryon as well as lepton number, and the process is found to be in 
equilibrium for a range of temperature i.e.,100Gey <T< lO^^GeV. 

For a systematic estimation of baryon asymmetry, one has to correlate it 
with the seesaw mechanism where both Dirac neutrino mass matrix and right- 
handed Majorana mass matrix are figured. In this context we refer to our 
earlier work[2] where we had employed the seesaw mechanism to successfully 
generate the degenerate, inverted and normal hierarchical models of neutrino 
mass patterns, using the diagonal texture of rriLn ^-nd non-diagonal texture 
of Mfiji. In most of the left-right symmetric SO(10)GUT models, the Dirac 
neutrino mass matrix mi^n can be either charged lepton mass matrix (referred 
to as case i) or up-quark mass matrix (referred to as case ii). However, we find 
that the seesaw mechanism alone fails to discriminate the above two possible 
choices of niLR, as well as the correct choice of the neutrino mass model in 
question. In the present work, we will estimate the baryon asymmetry of the 



universe using the same input values of the parameters which were aheady fixed 
for the predictions of masses and mixing angles at the seesaw stage. This in fact 
assures good predictions of neutrino masses and mixings. As emphasised before, 
such additional information on the calculations of the baryon asymmetry may 
be used to further discriminate the two choices of mLn, and also the correct 
pattern of neutrino mass model in question. 

In section 2, we present a brief review on the main formalism for estimating 
baryon asymmetry through the out-of-equilibrium decays of heavy right-handed 
neutrinos. Section 3 presents the numerical calculations of the baryon asym- 
metry for different neutrino mass models, and their comparison. Section 4 
concludes with a summary and discussion. Important expressions [2] related to 
rriLL and Man for all the neutrino mass models derived from seesaw formula, 
are relegated to Appendix A for ready reference in the present calculations. 

2 A brief review on the estimation of baryon 
asymmetry 

2.1 Type-I seesaw formula and lepton asymmetry 

Type-I seesaw formula[l] relates the left-handed Majorana neutrino mass matrix 
TO^^ and heavy right handed Majorana mass matrix Mjiji in a simple way 

■niLL = -miRMj^li^mlj^ (1) 

where tolr is the Dirac neutrino mass matrix. In some left-right symmetric the- 
ories such as SO(10)GUT, the left-handed Higgs triplet A^ acquires a vacuum 
expectation value and this leads two sources in the neutrino mass matrix[ll]: 
iTT'LL = "^LL + '^LL' where the first term is the usual seesaw term and the second 
one can be expressed as m£^ = "{{Mw/vRfMnR. In general 7 is a function of 
various couplings, and without fine tuning 7 is expected to be order of unity ( 
7 ~ 1) . The modified seesaw formula (referred to as Type-II) can be written 
as 

rriLL = -rriLRM^lmlR + "f{Mw/vRfMRR (2) 

Three different possibilities can arise between the two contributing terms i.e., 
to£^ and Ti^^. The situation where the second term is dominant over the 
first term, has the potential to connect [12,13,14,15] the large atmospheric neu- 
trino mixing and b-r unification in the context of minimal supersymmetric 
SO(10)models. In Type II seesaw formula, the lepton asymmetry ei can be 
written as: 

ei=e' + e" (3) 

where e^ and e^^ are the contributions of first term and second term of cq.(2) re- 
spectively. For our calculation of lepton asymmetry we consider the model[5,6,7] 
where the asymmetric decay of the lightest of the heavy right-handed Majorana 
neutrinos, is assumed. We confine our calculation to the Type I contribution 
only, assuming the extra contribution is relatively very sinall(7 ^ 10^*). 



ii:,M(Mt)y 

8- EJWhP ^^'' 
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The physical Majorana neutrino Nr decays into two modes: Nr -^ II + <P^ 
and II + 4>, where II is the lepton and II is the antilepton. For CP-violating 
decay through one-loop radiative correction by Higgs particle, the branching 
ratio for these two decay modes is likely to be different. The CP-asymmetry 
which is caused by the intereference of tree level with one-loop corrections for 
the decays of lightest of heavy right-handed Majorana neutrino A^i is defined 
by[6,8] e = ^ where T = r{Ni -^ IliJ)^) and f = r(iVi ^ hcf)) arc the decay 
rates. A perturbative calculation from the interference between tree level and 
vertex plus self-energy diagrams for A^i decay gives [6, 7, 16] 



ei 



h — mLji/v , and 

with j=2,3. The function fv{xj) represents the vertex contribution and fs{xj) 
represents the self-energy contribution of the decay process of the right-handed 
Majorana neutrino. The value of f{xj) in MSSM is twice its value for SM. For 
Xj >> 1, the expression (4) for SM case, can be written in terms of neutrino 
mass matrix as[17,18,19,20]: 

, _ _ 3Mi Im[{h*mi^h^U ... 

' - '^ - 167rt;2 {hh^),, ^^> 

For Xj ~ 1, we will have a very particular situation where Mi ~ M2 < M3, 
known as quasi-degenerate spectrum. The asymmetry is largely enhanced by a 
resonance factor which is given by[21,22,23]: 

^ |Mi| 

2(|M2|-|Mi|) ^°^ 

In such situation, we use the following expression [21] to calculate the lepton 
asymmetry, 

, ^ M2 Jm[(/i*mii/it)n] 



87rv2 (/i/it) 
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'-R (7) 



2.2 Prom lepton asymmetry to baryon asymmetry through 
sphaleron process 

The (B-L) conserving and (B+L) violating electroweak sphaleron processes [5, 6] 
give a set of coupled [24] equations for the evolution of lepton number and baryon 
number. At temperature T above the electroweak phase transition temperature 
(EWPT) T(7, the baryon asymmetry can be expressed in terms of (B - L) number 
density as [9]: 

24 + 4:Nm 



Again (B-L) asymmetry per unit entropy is just the nagative of til/ s^ since 
the baryon number is conserved in the right-handed Majorana neutrino decays. 
At the electroweak phase transition temperature, any primordial (B+L) will be 
washed out and the relation(8) can be written as[21,25]: 

nB_ ^ _24 + 4Ni^nL .. 

s ~ 66+13Nh s ^ ' 

where Nh is the number of Higgs doublet and til/s represents the lepton number 
to entropy ratio. For standard model(SM), Nh — 1; and eq.(9) reduces to 

^ c. -^!^ (10) 

s 79 s ^ ' 

For very slow B non-conserving interactions, the baryon to entropy ratio in a 
comoving volume is conserved. But the baryon to photon ratio does not remain 
constant due to the variation of photon density per comoving volume [8, 26] at the 
different epoch of the expanding universe. Considering the cosmic ray microwave 
background temperature T ~ 2.7 K, we have s — 7.04n^. For SM case we have 

S 

where ki is the dilution factor. This leads to the observed baryon asymmetry 
of the Universe [21], 

ySM = (^!^-^SM ^ ^^^^^ ^ 0.0216Kiei (11) 

For our numerical estimation of baryon asymmetry for different neutrino mass 
patterns, we use the above expression. 

2.3 Out-of-Equilibrium Decay 

To get the resultant baryon asymmetry, an out-of-equilibrium situation has to 
be provided by the expanding universe. The desired condition is satisfied if the 
temperature is smaller than the mass Mi of the decaying neutrino at the time 
of its decay. The inverse decay is blocked[6,8,9] if the expansion rate is greater 
than the interaction rate F. Thus, 

^ ^ 1.66VrT' 



Mpi 

where g* is the effective numbers of degrees of freedom available at the temper- 
ature T and Mpi is the Planck mass scale. Substituting the expression of F for 
the decay of Mi, we have, 

(fefet)^^Mi ^ IM^T^ 



For attaining out-of-equilibrium decay of Mi we equate T = Mi , and we use a 
parameter [27] derived from eq(12) as, 

(MtWl^ ^\^ (13) 

How much the produced asymmetry is washed out is described by Bohzniann 
equation and can be parametrized by another parameter k known as dilution 
factor[8,27]: 

, ^/ Kwk^ ifl0<i^<106 

"ni7#^ ifO<i.<10 (1^) 

An equivalent description for the out-of-equilibrium decay can also be expressed 
by the decay parameter 

K ~ mi/m* 

where the effective neutrino mass is defined as 

^ Ml 

and the equilibrium neutrino mass is given by [28]: 

Here, g* = 106.75, Mpi = 1.2 x lO^^GeV, v = 174GeV. Thus, for mi < m*, 
the heavy right-handed neutrino will satisfy the out-of-equilibrium condition. 
The values of K for different neutrino mass models under consideration, ranges 
from 3.7 x 10^ to 10~^ as shown in table-3. Therefore, the dilution factor ki 
which takes into account the washout process due to inverse decays and lepton 
number violating scattering, are calculated from expression(14). 

3 Numerical calculation and results 

To start with, we predict the left-handed Majorana neutrino mass ttlll using the 
seesaw forniula[2]. As emphasised earlier, we adopt the diagonal form of Dirac 
neutrino mass matrix ttilji for case(i) where the Dirac neutrino mass matrix 
is taken as charged lepton mass matrix and case(ii) where the Dirac neutrino 
mass matrix is taken as up-quark mass matrix. Using the non-diagonal form 
of right-handed Majorana mass matrix Mfm we obtain three patterns of neu- 
trino mass models: degenerate, inverted hierarchical and normal hierarchical. 
The corresponding expresions for ttill and Mfm for these models, are given in 
Appendix-A for ready reference to the present calculation. 

In the next step, the physical right-handed Majorana neutrino N' is defined 
in the basis where Mrr is diagonal with real positive eigen- values, i.e., 

N' = U]^N (16) 



where U]^ is extracted from the relation, 

Mrr' = UIMrrUr = Diag{Mi , M2 , M3) (17) 

This amounts to the re-definition of Dirac neutrino mass matrix ttilr- 

rriLR -^ m'j^R = tulrUr (18) 

In general Ur is complex, and the new m^^ in eq.(18) is made complex. It can 
be emphasised that such transformation makes the seesaw term toll invariant. 
In most calculations the Majorana CP phases enter through the construction 
of rriLL but in the present case such complex phases are originated from the 
diagonalising matrix Ur of Mrr. However both approaches are equivalent as 
they are related by an inverse seesaw relation Mrr = —niLRm'J^j^rrv^^^ which 
is specially true for diagonal form of mass matrix ttilr. In our calculation four 
neutrino mass models given in appendix viz., DegTlA, DegTlC, InvT2B, and 
NRT3, have complex diagonalising matrix Ur containing complex Majorana 
phases. But for two models DegTlB and InvT2A, the diagonalising matrix 
Ur is real, leading to zero lepton asymmetry. Since the non-zero lepton asym- 
metry requires complex Yukawa couplings i.e., complex phases, the following 
transformation is introduced by hand. Thus 

Ur ^ UrD (19) 

where, 

D = 

contains the complex Majorana phases. For further calculation of lepton asym- 
metry for the above two models we take a — P = tt/2. It can be emphasised 
here that such extra phase transformation of Ur does not change the seesaw 
term ttill and also the term {hh'') in cq.(4). 

For demonstration we first consider the normal hierarchical model as an 
example. The corresponding heavy right-handed neutrino mass matrix and 
light left-handed neutrino mass matrices are collected from ref. [2] where we use 
the diagonal form of Dirac neutrino mass matrix julr = i??ag.(A™, A", l)w. 
Normal Hierarchical Type3(NHT3): 

\ 2m— 1 \m-\-n — l \m — l 

Mrr = I A'"+"-l Xm+n-2 q j ^^ 

^m-1 Q I 



e*° 











1 











e^f' 



-A4 


A 


A--* 


A 


1-A 


-1 


A-^ 


-1 


1-A3 



V A-'' -1 1 - A3 y 

We take the input values A ~ 0.3, toq = 0.03ey leading to correct predictions 
of neutrino mass parameters and mixing angles: 



tan2 6*12 =0.55, sin^ 26I23 = 0.98, sin 613 = 0.074 

These are almost consistent with recent data[2]. The predictions on neutrino 
masses and mixing angles for all the models, are presented in Table- 1. 

For case(i): We use the charged lepton mass matrix, tue — itllr, {m,n) = 
(6, 2), vo = 1.01 X lO^s GeV, mo= 0.03 eV, tulr = tue = diag(X^, A^, l)i;,with 
V = 174GeV. This leads to 

1.7887 X 10" 2.2083 x lO" 2.4537 x lO^^ 
Mrr = I 2.2083 X 10" 7.3610 x 10" 

2.4537 X 10^2 1.0097 x lO^'"^ 

For real and positive eigen values [29] of Mfm the corresponding diagonalising 
matrix is 

0.964044e"*^ -0.265732 0.00242998 

Ur = I -0.265732e-'-5 -0.964047 5.31822 x lO"'^ 
-0.00234248e-^5 0.000646238 0.999997 

and the three heavy masses are 

M]^^^ = diag{6.51 x 10^°, 7.97 x 10", 1.01 x lO^^). 

The Yukawa coupling 'ft,' is defined[23] hy h — mLu/v, where v is the electroweak 
vacuum expectation value. As the Dirac neutrino mass matrix is rotated in the 
basis where right-handed Majorana neutrino mass matrix is diagonal, the new 
Yukawa coupling becomes: h — mLjiUu/v 

-0.000702788e-*f -0.000193718 1.77146 x 10"^ 
0.0239159e-'t -0.0867642 4.7864 x 10"^ 

0.00234248e-*^ 0.000646238 0.999997 

For case(ii): We again use the up-quark mass matrix niup = rn-LR, {m,n) = 
(8,4), Vo = 1.01 X lO^^GeV, too = 0.03 eV , tulr = TO„p = diag{X^,X^,l)v. 
This leads to 



1.4495 X 10^ 1.7895 x 10^ 2.2092 x 10" 

RR 



Mrr = I 1.7895 X 10^ 5.9648 x 10^ 



2.2092 X 10" 1.0102 x lO^'"^ 

The corresponding diagonalising matrix is: 

0.964047e-*f -0.265733 0.0002187 

Ur = I -0.265733e-'t -0.964047 -3.87423 x 10" 

-0.000210837e-^* 0.0000581162 -1 
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Type 


Am2jlO-5ey2] 


Am23[10-3ey2] 


tan^ 6i2 


sin^ 2023 


sin 013 


DegTlA 


8.80 


2.83 


0.98 


1.0 


0.0 


DegTlB 


7.91 


2.50 


0.27 


1.0 


0.0 


DegTlC 


7.91 


2.50 


0.27 


1.0 


0.0 


InvT2A 


8.36 


2.50 


0.44 


1.0 


0.0 


InvT2B 


9.30 


2.50 


0.98 


1.0 


0.0 


NHT3 


9.04 


3.01 


0.55 


0.98 


0.074 



Table 1 : Predicted values of the solar and atmospheric neutrino mass-squared differences and 
three mixing parameters calculated from »Tt£^ derived from seesaw formula in the Appendix 
A. 



M^^^ = diag{5.2699 x 10^ 6.4571 x 10^ 1.01 x lO^^) 



The Yukawa coupling 'h' is 



h = 



-7.02787 X 10-''e-*f 


-1.93719 X lO-'' 


1.59431 X 10-^ 


0.02391e-'5 


-0.08676 


3.4868 X 10-1^ 


2.10837 X lO-'^e-** 


5.81161 X 10-5 


1 



To calculate the dilution factor first we calculate the effective mass m for each 
case (Table-3). For normal hierarchical case we have for case(i): rfii = 2.47 x 
lO""* eV and for case(ii): mi = 2.47 x 10^^ cV The corresponding dilution 
factors are 0.17 and 0.17 for case(i) and case(ii) respectively. Considering the 
decay of A'^i, the baryon asymmetry is found to be: 
for case(i) where itilr = rnE, 

Yi^^ ^ 2.17 X 10"^ 

and for case(ii) where rriLR — niupjY^^'' — 1.76 x 10~^^ Following the same pro- 
cedure we calculate the baryon asymmetry for all the models given in Appendix- 
A. The results are presented in Tables 1-5. As emphasised earlier, we use only 
the values of the parameters fixed at the first stage where seesaw formula oper- 
ates, in the estimation of baryon asymmetry in the present calculation. 

4 Summary and Discussion 

We summarise the main points in this work. We have attempted to give a comparison 
of the numerical estimations on baryon asymmetry of the universe for different neutrino 
mass models. We have investigated these neutrino mass models (degenerate, inverted 
hierarchical and normal hierarchical) at two consecutive stages: (i) from neutrino 
oscillation parameters (Table 1) using seesaw mechanism and (ii) from the observed 
baryon asymmetry (Table 4), within a unified framework. In a model where baryon 
asymmetry is explained via lepton asymmetry through sphaleron process, the out-of- 
equilibrium decay of the right-handed Majorana neutrino plays a decesive role. As the 
out-of-equilibrium decay of the lightest heavy Majorana neutrino Mi is essential for 



Type 


Case(i): Mj 


Case(ii):|Mj| 


DcgTlA 


4.28 X 10■^ 1.16 X 10^", 3.84 x 10^^ 


3.47 X 10^ 9.42 x 10', 3.81 x lO^"^ 


DegTlB 


4.05 X 10^6.16 X 10", 7.6 x 10^^ 


3.28 X 10^ 4.98 x 10^ 7.6 x 10^^ 


DegTlC 


4.05 X 10^6.69x 10i2,6.99x lO^^ 


3.28 X 10^,4.85 X 10", 7.81 x 10" 


InvT2A 


3.28 X 10^9.70 X 10^^^6.79 x lO^^ 


2.64 X 10^, 7.92 x 10^°, 6.70 x lO^^ 


InvT2B 


5.6527 X IQi", 5.6532 x 10^°, 5.38 x lO^^ 


4.5971 X 10^,4.5974 x 10^,5.34 x lO^^ 


NHT3 


6.51 X 10^", 7.97 X lO^i, 1.01 x lO^^ 


5.27 X 10^6.45x lO^l.Olx lO^'"^ 



Table 2: The three right-handed Majorana neutrino masses in GeV for both case (i) and 
case (ii). The expressions of Mnn are collected from Appendix A. 





For 


casc(i) 




For 


case(ii) 




Type 


rhi 


K 


Kl 


mi 


K 


Ki 


DcgTlA 


3.76 X 10-^ 


3.48 


0.11 


3.74 X 10-3 


3.46 


0.11 


DegTlB 


0.40 


370 


2.79 X lO-'^ 


0.40 


370 


2.79 X lO-'^ 


DegTlC 


0.40 


370 


2.79 X lO-'^ 


0.40 


370 


2.79 X lO-'^ 


InvT2A 


0.05 


46.3 


2.9 X 10-3 


0.05 


46.3 


2.9 X 10-3 


InvT2B 


2.85 X lO-'^ 


0.26 


0.17 


2.83 X 10-* 


0.26 


0.17 


NHT3 


2.47 X lO"'^ 


0.23 


0.17 


2.47 X 10-* 


0.23 


0.17 



Table 3: Calculation of values of effective mass parameter mi in eV, and dilution factor ki. 



Type 


e^ (case i) 


€^ (case ii) 


Y^'^'Caseii) 


Y^'^'Caseiii) 


DegTlA 


2.10 X 10-" 


1.71 X 10-» 


4.99 X 10-^ 


4.06 X 10-" 


DegTlB 


2.66 X 10-1^ 


2.16 X 10-2" 


1.60 X 10-23 


1.30 X 10-25 


DegTlC 


1.74 X 10-**^ 


1.69 X 10-2" 


1.05 X 10-23 


1.02 X 10-25 


InvT2A 


1.59 X 10-" 


1.27 X 10-1" 


9.94 X 10-19 


7.96 X 10-21 


InvT2B 


1.47 X 10-2 


1.62 X lO-'^ 


5.40 X 10-'"^ 


5.94 X 10-^ 


NRT3 


5.90 X lO-'^ 


4.78 X lO-'^ 


2.17 X lO-'^ 


1.76 X 10-" 



Table 4: Calculation of Icpton asymmetry e^ and baryon asymmetry Yg for casc(i) and 
case(ii) for respective neutrino mass models given in Appcndix-A. 



(to, n) 



(8,4) 



5.27 X 10« 
4.78 X 10-9 
1.76 X 10-" 



(6,2) 



6.51 X 10*" 
5.9 X lO-'^ 
2.17 X 10-9 



(4,2) 



1.16 X 10^2 

1 X 10-5 
3.24 X IQ-s 



(7,2) 



1.87 X 10*" 

1.7 X 10-^ 

6.24 X 10-1° 



Mi{GeV) 
ei 

VSM 
^ B 



Table 5: Calculation of baryon asymmetry for normal hierarchical mass model for different 
choices of (m, n) pair. 



fruitful leptogenesis, we consider the effective mass parameter mi for all the models 
presented in Table 3 where only inverted hierarchical model(InvT2B) and normal 
hierarchical model (NRT3) satisfy the out-of-equilibrium decay condition. In these 
two models the effective neutrino mass is less than the equilibrium neutrino mass 
m*(~ 10"'^) i.e., mi < m*. 

We present the numerical predictions on solar and atmospheric neutrino mass- 
squared differences and three mixing angles, which are calculated using seesaw formula 
in Table 1. All the input parameters in Mrr and niLR are fixed at this stage. Table 
1 indicates good predictions for normal hierarchical model and inverted hierarchical 
model( InvT2A). However, the inverted hierarchical (Type 2A) model is not stable 
under radiative corrections in MSSM unlike normal hierarchical model(NRT3) . In Ta- 
ble 2 we present the masses of the physical right-handed Majorana neutrinos, which 
are extracted through the diagonalisation of Mb.r- Further, the calculated values of 
baryon asymmetry for all neutrino mass models are shown in Table 4 which clearly 
shows the competitive nature of degenerate (Type lA) model [30] and Normal hier- 
archical model. For both models, the produced baryon asymmetry is slightly lesser 
than the experimental value in case (ii) and slightly larger in case (i). This possibly 
implies that the actual form of Dirac neutrino mass matrix may lie between these 
two extreme cases. A result of our investigations in this regards for normal hierar- 
chical mass model, is presented in Table-5 for different values of (m, n) pair in rriLR 
[31]. As the mass of the lightest right-handed Majorana neutrino Mi is concerned, 
the model agrees with the famous Ibarra-Davidson bound [32] i.e., Mi > 4 x 10* GeV. 
Considering the stability criteria under radiative corrections in MSSM as well as good 
predictions on oscillation mass paramerters and mixing angles, the normal hierarchi- 
cal model is quite favourable. Future neutrino oscillation experiments may confirm 
this conjecture. Finally we also point out some ambiguities on the choice of the value 
of coefficent d = 0.216. The present value is found to be 2.2 times larger than the 
value appeared in refs.[16,19,23,33] and 0.74 times smaller than the value appeared in 
ref.[34]. 

Appendix A 

Here we collect the various right-handed Majorana mass matrices Mrr and the light 
left-handed neutrino mass matrix rriLL from ref. [2] for ready reference in the text. 
Here we use diagonal form of Dirac neutrino mass matrix rriLR = Diag{\"^, A", l)v. 
The nomenclature of the different neutrino mass models has been changed here. 

Degenerate TypelA(DegTlA) with mass eigenvalues (mi, — m2,m,3): 




V2 



-5i)A'"+" {l/2 + Si^52)\^" {~l/2 + Si-S2)X" 



vo 



(^+5i)A'" (-l/2 + 5i-52)A" (-1/2 + 51-^2) 



(-251+252) (^-<5i) (75 --^i) 
-mL=( (75-^1) {1/2 + S2) (-1/2 + 52) I mo 

(^-5i) (-1/2 + 52) (1/2 + 52) 
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Here, 5i = 0.0061875, S2 = 0.0031625, A = 0.3, m,, = 0.4eV and vq = 7.57 x 10"GeV 
Degenerate TypelB(DegTlB) with mass eigenvalues {mi,m2,m3): 

( (1 + 25i + 252)A2'" (5iA"+" (5iA" \ 

Mrr= 5iA™+" (l + ^aA^" &.,r- K;o 

V (5iA™ 52A" (l+<52) / 

/ (1 - 2<5i - 2^2) -&i -5i \ 

-niLL = -Si (l-<52) -<52 mo 

V -<5i ~S2 {I-S2) J 

Here, 5i = 7.2 x 10"^ <52 = 3.9 x 10"^ A = 0.3, rm, = 0.4cV and vo = 7.57 x lO^GeV. 
Degenerate TypelC(DegTlC) with mass eigenvalues (mi, 7712, —ms): 

/ (1 + 25i + 2(52)A2'" <5iA'"+" 5iA™ \ 

Mrr= 5iA'"+" 52A2" (1 + 52)A" Lo 

V ^lA'" (1 + <52)A" 52 J 

/ (1 - 25i - 2^2) -5i -5j \ 

-hiLl = -5i -(52 (1 - <52) mo 

V -5i (l-'52) -<52 / 

Here, 5i = 7.2 x 10"^ (52 = 3.9 x 10"^ A = 0.3, mo = 0.4cV and «o = 7.57 x lO^^GeV. 
Inverted Hierarchical Type2A(lnvT2A) with mass eigenvalues (m,i, m2,m,3): 



Mrr= rie\"'+" 1/2A2" -(1/2 - rj)\" \^ 

\ rjeX"' -(l/2-?7)A" 1/2 / '' 



(l-2e) -e -e 

i (- 

2 V2 



-m^L = -e 5 (I -»?) 1 mo 



Here, »? = 0.0045, e = 0.0055, A = 0.3, mo = 0.05eV and vo = 6.06 x lO^GeV 
Inverted Hierarchical Type2B(lnvT2B) with mass eigenvalues (mi,— m,2,m3) 

C\2m + 7 \Tn+n+4 \m+4 \ 

^m+n+4 ^2n _^n \ ^ 

^m+4 _-j^n -^ / 
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-A" 


A 


A^ 


A 


1-A 


-1 


A^ 


-1 


1-A=* 



/ 1 1 \ 

-mLL= 1 -(A3-A*)/2 -(A« + A4)/2 mo 
V 1 -(A3 + A^)/2 -(A3-A4)/2 / 

Here, A = 0.3, mo = 0.035eV and wo = 5.34 x 10^'^GeV 

Normal Hierarchical Type3(NHT3) with mass eigenvalues (7711,7712,7713) 

(\2m — 1 \7n+n — 1 \m— 1 \ 

^m + n^l ^m + n-2 q j „^ 

A™-' 1 / 



-hiLl =1 A i - A -i I mo 



Here.A = 0.3, mo = 0.03eV and ^o = 1.01 x lO^'^GeV. 
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